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ABSTRACT 

Redshifted HI 21 cm emission from unresolved low-redshift large scale structure is a promis¬ 
ing window for ground-based Baryon Acoustic Oscillations (BAO) observations. A major 
challenge for this method is separating the cosmic signal from the foregrounds of Galactic 
and extra-Galactic origins that are stronger by many orders of magnitude than the former. The 
smooth frequency spectrum expected for the foregrounds would nominally contaminate only 
very small A:|| modes; however the chromatic response of the telescope antenna pattern at this 
wavelength to the foreground introduces non-smooth structure, pervasively contaminating the 
cosmic signal over the physical scales of our interest. Such contamination defines a wedged 
volume in Fourier space around the transverse modes that is inaccessible for the cosmic sig¬ 
nal. In this paper, we test the effect of this contaminated wedge on the future 21 cm BAO 
surveys using Fisher information matrix calculation. We include the signal improvement due 
to the BAO reconstruction technique that has been used for galaxy surveys and test the effect 
of this wedge on the BAO reconstruction as a function of signal to noises and incorporate the 
results in the Fisher matrix calculation. We find that the wedge effect expected at z = 1 — 2 
is very detrimental to the angular diameter distances: the errors on angular diameter distances 
increased by 3-4.4 times, while the errors on H{z) increased by a factor of 1.5-1.6. We con¬ 
clude that calibration techniques that clean out the foreground “wedge” would be extremely 
valuable for constraining angular diameter distances from intensity-mapping 21 cm surveys. 

Key words: cosmology: observations, 21cm, distance scale, large-scale structure. 


1 INTRODUCTION 


Baryon acoustic oscillations (BAO) imprinted in the large scale 
clustering of matter have emerged as one of the leading methods to 
measure the effect of the mysterious dark energy on the geometry 
of the Universe (e.g. Hu&Whhe|l^% Eisenrieinj200^ Blake & 


Glazehrook|2003l|Linder|2003||Hu & I aiman |2003||Seo & Risen 

stein|2003) . BAO were formed hy primordial sound waves that had 
propagated through tightly coupled photons and haryons in very 
early Universe and subsequently freezed out at the epoch of re¬ 
combination when the photons and haryons decouped. The distance 
that the sound wave has traveled before the recombination epoch is 
called the sound horizon scale, and this sets the physical size of this 
feature. The sound horizon scale straightforwardly depends on the 
sound speed and the time the sound wave traveled, which are pre¬ 
cisely calibrated by the cosmic microwave background data. There¬ 
fore the observed sizes of BAO in comparison to this true physical 
size construct a robust standard ruler test, providing cosmological 
distances, such as angular diameter distances and Hubble parame¬ 
ters to various redshifts, i.e., the expansion history of the Universe. 

The sound horizon scale today is ~ 150 Mpc and BAO sur¬ 
veys have been designed to measure the large scale structure over 
many hundreds of Mpc or a few Gpc in order to reduce statis¬ 
tical sample variance in the measurements of the 150 Mpc-long 


BAO scale. BAO surveys historically have focused on the opti¬ 
cal band using spectroscopic or photometric methods, where we 
measure the large scale clustering in galaxies (or quasars and LyA 


forests) as biased tracers of the underlying dark matter (e.gJColc 


et al.|2005l|Eisenstein et al.|2005l|Percival et al.|20071|2010l|Blake 

et al.|201 la|bl IBeutler et al.||201 11 iPadmanabhan et al.||2012l |Xu 

et al.||2012l ISeo et al.||2012l IKazin et al. 

120131 lAnderson et al. 

20141 IKazin et al.|2014l|Font-Ribera et al. 

2014l|Ross et al.|20151 

Delubac et al.|2015^. 

tral hydrogen has long 
'reionization. However, 
e found in galaxies, and 
jresents a possible route 
f the BAO signal (e.g.. 

Redshifted 21 cm emission from neu 
been investigated as a probe for the epoch o 
even after reionization neutral gas can still h 
thus intensity mapping in the 21 cm line re] 
to low and intermediate-z measurement o 

Peterson et al.|2006 Chang et al.|2008||Ansari et al.|2008||Wolz 


et al.||2014[ [Bull et al.||2015|. Among many potential advantages 


of the 21cm surveys, we note the large instantaneous field of view 
and hence fast survey speed, comparatively low cost of the tele¬ 
scope, feeds, and electronics, a spectral (hence redshift) resolution 
that can be very high since it is determined by digital signal pro¬ 
cessing rather than optics, and absence of the complex atmospheric 
absorption and emission features that present a challenge in the far- 
red spectral domain. Such advantages and promising scientific re¬ 
turn have motivated ongoing efforts for near future 21cm BAO ex- 
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2 Seo & Hirata 


periments such as the Canadian Hydrogen Intensity Mapping Ex¬ 
periment (CHIME) [5 the Baryon Acoustic Oscillation Broadband 
and Broad-beam (BAOBAB, |Pober et al.|2013^ , BAORadio ( |Ansari| 
|et al.|[20T2) , the BAO from Integrated Neutral Gas Observations 
experiment (BINGO, IBattye et al.|20I3), Ooty dAIi & Bharadwajl 
|20T4| > and others. 

The biggest challenge of this method is separating the as- 
trophysical foreground emissions from the cosmic signal that is 
weak by many orders of magnitudes than the former. The good 
news is that the foreground due to Galactic and extragalactic syn¬ 
chrotron and free-free emission are spectrally smooth, which means 
they in principle would only contaminate very small line-of-sight 
wave modes in Fourier space (i.e., mainly fluctuations in the trans¬ 
verse direction). Losing the small line-of-si ght Fourier modes (i.e., 
fcii < 0.01 — 0.02h~^ Mpc, references in Chang et al. 


2008 


has 


very little impact on the BAO error forecasts and therefore many, 
if not all, error forecasts for the 21cm BAO surveys assume a per¬ 
fect foreground subtraction taking this advantage and include the 
foreground effect mainly in the error budget (e.g.,|Seo et al.|2010bl 
[Bull et al.|20T5l ). 


In reality, however, the frequency dependence of an interfer¬ 
ometer’s fringe pattern (i.e., chromatic instrument response) can 
cause the spectrally smooth foreground spectrum to ring at a dif¬ 
ferent angular location as a function of frequency. This introduces 
an un-smooth foreground component and therefore makes a clean 
separation between the foreground and the cosmic signal very dif¬ 
ficult (e.g., [Alonso et al.||20l3) . However, it does turn out to be 
confined to a so-called “foreground wedge” at low fc|| and high k± 
( jPatta et al.||2010l jVedantham et al.||2012[ [Morales et al.|[20T2l >, 

leaving a cleaner “window” at high A:|| and low k±. The wedge ef¬ 
fect is expected to be much alleviated for surveys with very small 
side lobes, such as single dish surveys with low aberrations and un¬ 
obstructed apertures (e.g. the surveys performed at the Green Bank 
Telescope ( [Masui et al.[2013[ l). For interferometer arrays, side lobes 
due to various instrumental effects that scale with frequency will 
introduce the foreground wedge. The wedge is not a fundamental 
effect, in the sense that with excellent baseline-to-baseline calibra¬ 
tion the foregrounds in this region can be removed, but the factor by 
which raw systematics dominate over cosmological signal is orders 
of magnitude greater in the wedge than the window (see e.g. |Th^ 
[garajan et al.|2015l for a discussion in the context of epoch of reion¬ 
ization studies). It is thus of great interest to understand the impact 
of the wedge on cosmological constraints. How are cosmological 
constraints affected if we assume that only the cosmological signal 
in the foreground “window” is available? Or, stated more opti¬ 
mistically, what are the potential gains to cosmology if advances 
in side-lobe mapping and data processing techniques can clean out 
the wedge? [Pober] ( [2015^ , for example, has shown the impact of 
the foreground wedge on redshift-space distortion measurements, 
which is destructive at a ~ 8 but can place non-trivial constraints 
at 2 ~ 1. 


The purpose of this paper is to estimate the effect of remov¬ 
ing the wedge on the BAO analysis. We use the Fisher information 
matrix calculation as well as A-body simulations to consider the 
wedge effect on the BAO reconstruction. 

The expected performance of a 21 cm BAO survey in the pres¬ 
ence of the wedge effect could be straightforwardly calculated us¬ 
ing the Fisher matrix formalism by removing the contribution from 
the Fourier modes in the wedged volume. Such calculation has a 


^ http://www.mcgillcosmology.ca/chime 


very simple dependence on the expected sample variance and the 
expected strength of the BAO signal from the assumed large-scale 
clustering data. The strength of the BAO signal is reduced at low 
redshift due to various nonlinearities in the process of structure 
growth and redshift-space distortions in a predictable manner. For¬ 
tunately, a large portion of such nonlinear reduction of the signal 
can be undone by using the BAO reconstruction technique ( [Eisen-[ 
[stein et al.[2007b^ . In this technique one uses the observed cluster¬ 
ing of a given tracer to estimate the gravitational potential field due 
to the underlying matter density field and use the estimated poten¬ 
tial field to correct for the motions of mass elements that weakened 
the BAO signal. The BAO reconstruction is proven to be successful 
in subsequent simulations (e.g., Huff^^dJ2007 Seo^^e^^an2008 


Padmanabhan et al.[2009||Noh et al.[2009|[Seo et k.[2010a[|Mehta 


et al.|20lT|l as well as in recent galaxy redshift surveys (e.g., [Pad^ 


manabhan et al.[2012[ [Anderson et al.[2012[>, and now has become 


a standard tool for BAO surveys with a good redshift precision. 

In order to include the effect of the BAO reconstruction in the 
Fisher matrix formalism, we need an estimate of the strength of 
the BAO signal after reconstruction. In a large signal-to-noise limit 
(per Fourier mode or per pixel), we conventionally assume that the 
width of the BAO peak reduces to about a half after reconstruction 
at 2 ~ 0.3 ( [Eisenstein et al.[2007b[|Seo & Eisenstein[2007[ l. 21cm 
BAO surveys, however, operate often at a low S/N per mode due 
to the high system (foreground as well as instrument) temperature 
relative to the weak line signal, while the total S/N contributed by 
many modes is boosted by the large survey area. The reconstruction 
efficiency is in general expected to decrease as the noise level of the 
density field increases; [Seo et ^ ( [2008^ and [Mehta et al.|j2011[ l 
used simulations to show that a density field with S/N per mode 
(at k = 0.2h Mpc“^) as low as ~ 0.6 nevertheless shows a non- 
negligible improvement by the reconstruction; recently, [Burden ^ 
[ar] ( [2014| l tested a regime with S/N > 0.2 for galaxy surveys and 
found a similar result. In the current paper, we test a density field 
with S/N as low as 0.1. 

An analytic prediction on the BAO reconstruction perfor¬ 
mance as a function of (shot) noise level has been derived based 
on the Lagrangian Perturbation Theory formalism by [White] j2010[ > 
in real space. The presence of the wedge contamination makes esti¬ 
mating the effect of the BAO reconstruction using such an analytic 
approach nontrivial especially when including redshift-space dis¬ 
tortions. In this paper, we investigate the effect of the wedge on the 
BAO reconstruction using A-body simulations as a function of a 
signal-to-noise ratio per mode as low as ~ 0.1. We incorporate our 
results in the Fisher matrix calculation and derive the final effect 
of the wedge contamination on the distance precisions expected for 
21 cm BAO surveys. 

This papers is structured as follows. In § |^ we set up Fisher 
matrix calculations in the presence of the wedge effect as well as 
introducing a fiducial 21cm BAO survey that we adopt from [Seo] 
[et al.[ j2010b[ l. In § [^ we present the details of performing BAO 
reconstruction in the presence of the wedge at different redshifts 
and at various signal to noise levels that are simulated. In §|^ we 
present the results of this paper, the BAO reconstruction efficiency 
at different signal-to-noise levels as well as at different extents of 
the wedge, and the Fisher matrix calculations incorporating BAO 
reconstruction. Finally, in §|^ we conclude. 


© 2011 RAS, MNRAS 000,[T]-?? 








































































The foreground wedge and 21 cm BAO surveys 3 


2 SETTING UP THE EISHER MATRIX ANALYSIS EOR 
21CM BAO SURVEYS 

The Fisher matrix formalism has been an effective tool for pre¬ 
dicting constraints on cosmological parameters for future BAO sur¬ 
veys. In this section, we lay out the components of the Fisher matrix 
calculations aimed at incorporating the wedge effect. We start from 
the method in |Seo & Eisenstem| ( |2007| l that isolates the BAO-only 
information. In order to quantify the effect of the wedge, we adopt 
the same fiducial 21cm survey configuration described in |Seo et 
( |2010b^ and compare the results before and after incorporating the 
wedge effect. 

2.1 Fiducial cosmology of the Fisher analysis 

For simplicity, we define the fiducial cosmological parameters to 
be = 0.27, h = 0.72, Oa = 0.73, Qk = 0, = 0.0238, 

T = 0.17, Us = 0.99, which is one of the cosmologies uset^jin 
|Seo & Eisenstein| ( |2007[ >, rather than using the best fit cosmologies 
from the most up-to-date CMB probes. We are therefore testing the 
effect of the wedge given this cosmological assumption; however, 
we do not believe that the relative effect of the wedge will sensi¬ 
tively depend on this assumption. For example, the wedge angle is 
not much different between the fiducial cosmology and the Planck 
{TT spectra -l- polarization at low multipoles) best fit ACDM cos¬ 
mology ( [Planck Collaboration et al.|2015^ ; the cosine of the angle 
is different only by 3% (Figure 1). 


2.2 Fiducial 2Icm survey configuration 

When investigating the wedge effect, we assume a fiducial 21cm 
survey in order to approximately simulate a reasonable redshift 
range as well as the signal-to-noise range to consider. We adopt 
the fiducial CRT survey assumed in jSeo et aL| ( |2010b^ that rep¬ 
resented one possible example of cost-effective future 21cm BAO 
survey designs. Table[^in this paper is taken from Table 1 in|Seo| 
jet al.|p010b [( that defines the fiducial survey configuration. While 
we refer interested readers to jSeo etn^ ( |2010b) for more details on 
assumptions, we summarize the key aspects here. 

The signal to noise of the BAO feature in the 21cm inten¬ 
sity power spectrum when averaged over dk and d^, where fc is a 
wavenumber and is the cosine of the line-of-sight angle, depends 
on the strength of the 21cm signal due to the large scale clustering 
of the neutral hydrogen as signal and instrumental noise and the 
sky background!^ 


s 

/ 27 rfc 2 dk dfj. 14urvey 

PHi(fc, fJ.) lU^ 

A ^ 

/ 2(277)3 

Phi(^, + Pth{k, fl) + Pshot 


( 1 ) 


where Kurvey is the total volume of the survey and we add a galac¬ 
tic shot noise contribution Pshot = l/h- due to the discreteness 
effect of the HI sources with an effective number density n ^Seoj 
jet al.|2010b^ . Sxy and T,z are the Gaussian damping lengths that 


^ This was based on the first year Wilkinson Microwave Anisotropy Probe 
jSpergel et al.|2003| 

^ We include here the BAO damping factor that is not explicit in the equa¬ 
tion in|Seo et al.|j2010b|. 

^ we assume a Gaussian density field; Takahashi et al. (2011) and Ngan et 
al. (2012) have shown that the non-Gaussian en'ors on the BAO measure¬ 
ment in multi-parameter fitting are negligible. 


Table 1. Fiducial CRT configuration from Table 1 of |Seo et ^j2010b) . For 
technical reasons, it was assumed that the redshift range is covered by two 
distinct configurations of the telescope. The fiducial telescope is a square in 
the sense that the total width -/Vcyi Wcyi is equal to Lcyi- 


Parameter Low redshift High redshift 


Redshift range 
Number of cylinders, A^cyi 
Length of Cylinder, Lcyi (m) 
Feed spacing, dp (m) 

Width of Cylinder, VFcyi (in) 
Duty factor, D f 
Wyear (years) 

bias 

Sky temperature, Tgky (K) 
Antenna temperature, Ta (K) 
gain, g 

^shot 


0.66 - 1.24 

1.22 - 2.11 

7 

10 

99.8 

142.8 

0.39 

0.558 

14.3 

14.3 

0.5 

0.5 

1.40 

0.87 

0.00037 

0.00037 

1.0 

1.0 

10 

10 

50 

50 

0.8 

0.8 

100.0 

100.0 


account for the degradation of the BAO signal strength due to var¬ 
ious nonlinearities. Due to redshift-space distortions along the line 
of sight, T>z tends to be larger than Exy. The values of Exy and Ez 
can be fairly precisely predicted as a function of cosmology and 
redshift in linear theory (e.g., |Eisenstein et al.||2Q07a[ |Matsubara| 
[2008| t. 

For a dual-polarization cylinder telescope built near the Equa¬ 
tor and with the long axis of the cylinders oriented North-South, the 
instrument thermal noise term, expressed in conventional galaxy 
power spectrum units (h~^ Mpc^), is 


Pthik,fL) 


/ ffTsky + Tg \ ^ dV 27rAobslT'cyl 
\ yTsig / dAdf year 


where Tgig is the average brightness temperature due to the 21cm 
line, Tsky is the average sky temperature (e.g., due to foregrounds), 
Ta is the average antenna noise temperature or the amplifier noise 
temperature, and Aobs = c(l -I- 2)//2icm is the observed wave¬ 
length jSeo et al.|2010b i. The conversion factor dV/dA df is from 
observed angle-frequency volume element to comoving cosmic 
volume, and has units of {h~^ Mpc)®/sr/Hz. This accounts for 
the distribution of time over the full range of right ascension (27r) 
and the declination range given by the Nyquist limit of the beam, 
I sinDecI < Aobs/ i^dp), where dp is the feed spacing; the survey 
area is 


Asurvey = 27r ^ = 271 Aobs , 

ap -Lcyl 


(3) 


so long as dp > Aobs/2 (as is the case here - otherwise the tele¬ 
scope is sensitive all the way to the celestial poles). The number 
of feeds is Afeed = Lcyi/dp. We have incorporated a duty factor 
D/ ^ 1 (the fraction of time the telescope is operating in science 
mode) and a survey duration Ayear. 

We assume that the neutral hydrogen traces the matter, so that 
Piiiik, fi) = 6^(1 + j3f^fPm{k) with bias 6=1 and redshift- 
space distortion parameter (Kaiser 1987 1 j3 = d In G(a)/ci In a/6 
where G(a) is the growth function. Note that assuming 6 > 1 im¬ 
proves the signal to noise. The signal temperature is (IBarkana &| 
|Loeb|2007[ [Pritchard & Loeb|2008[|aiang et al.|2008|l 


- 188 h{z)/Ho 


(4) 
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4 Seo & Hirata 


where xm is the neutral hydrogen fraction at z, SIh.o is the ratio 
of the hydrogen mass density to the critical density at 2 = 0, and 
H (z) is the Hubble parameter at 2 . 

The signal is suppressed by an instrument response functiorQ 
which will be modified in the next section to include the wedge 
effect. The general response function is 


W{k) 


^(wAobS) uAobs) 


sin(TA//2) 

rAf/2 


(5) 


where (w, v, r) is the Fourier mode in observed (angle-frequency) 
space that corresponds to k (the conversions are given explicitly in 
& and A/ is the width of a frequency channel (i.e. the pixel 
transfer function in the radial direction). The first term is the dis¬ 
tribution of baselines at separation (uAoba, vAobs), normalized to 
A(0, 0) = 1. We assume that the frequency channels are narrow 
(<C 1 MHz) so that the last term can be neglected. 

This response function is slightly different from the form 
adopted in |Seo et aL|p010bt . |Seo et al.|p010b^ assumed a specific 
choice of pixelization along the declination direction and accounted 
for the effect of the default number of cylinders that fill the baseline 
along the right ascension; here we assume a generalized response 
function approximation in both directions for simplicity. This gives 


Winat(fc) = (1 - |fcx/fcmax|)(l “ |A:y/fcmax|) (6) 


for fcx, ky ^ fcniax, and otherwise 0. Here fcx is a wavevector com¬ 
ponent in the direction of right ascension, ky is a component in the 
direction of declination, and fcmax is set to be 2'KLcyi/Xoba/D{z), 
where D{z) is the comoving angular diameter distance. For exam¬ 
ple, fcmax = 0.35(i Mpc“^ at 2 = 1 and 0.6h Mpc“^ at 2 = 2 
at our fiducial cosmology. Although not presented in this paper, we 
repeated the entire analysis in this paper with the response function 
presented in |Seo et al.| ( |2010b) ; qualitatively we find very similar 
results, while quantitatively the current response function produces 
larger errors[^As a caveat, although we do not put a limit on fcmax, 
all our Fisher matrix calculations do not account for a BAO infor¬ 
mation beyond 0.5h Mpc“^. 


2.3 The wedge effect due to the foreground contamination 

The wedge effect becomes important when we lift the assumption 
that the side lobes are negligible. We first discuss some considera¬ 
tions of the foreground before the assumption is lifted. 


The foreground subtraction for all cosmological 21 cm exper¬ 
iments - whether low -2 for BAO or high -2 for the epoch of reion¬ 
ization - depends on the spectral smoothness of synchrotron and 


free-free emission l|Di Matteo et al.|2004 

Gnedin & Shaver 

2004 

Morales & Hewitt 2004||Zaldarriaga et al. 

2004| Santos et al. 

2005 

McQuinn et al. 12006 1 , as compared to the cosmological signal for 


which the large-scale clustering is expected to introduce variations 
in the HI density and temperature as a function of observed fre¬ 
quency (or distance from the observer). Foreground mechanisms 


® This response function is often included in the noise rather than the sig¬ 
nal, in which case the noise power is inversely proportional to the number 
of baselines at the relevant separation. Where one should put it depends on 
whether the sky map under consideration is instrument-convolved or not; 
the choice has no impact on the signal-to-noise ratio or cosmological pa- 
ra meter constraints, so long as it appears somewhere! 

® |Seo et al.|(20T0bt adopted 

Winst(fc) = exp [-1.5(fcxAmax)^] ©(^max “ ky) 
with fc^ax = ^f-cyl/Aobs/TiCz)- 


with strong frequency-dependent structure, such as extragalactic ra¬ 
dio recombination lines, are expected to be small (e.g. |Oh & Mack| 
|2003| l. This phenomenon can be quantified by considering the cos¬ 
mological data cube in comoving coordinates {r^, ry,rz) (with Xz 
being the line-of-sight or redshift direction) and mapping it to the 
directly observable quantities {0x,0y,f) (with 9 being the angu¬ 
lar direction on the sky in orthographic projection with the primary 
beam line of sight at 6^ = 9y = 0, and / as the observed fre¬ 
quency; / = / 21 cm/(1 + z), where / 21 cm = 1-42 GHz is the rest- 
frame frequency of the 21 cm line). Therefore there is a Jacobian 
mapping between these spaces, which is locally diagonal: 


OTz; 

dOz, 


dry 

dOy 


D{z) 


(7) 


and 

dVz _ -c(l-|-2)^ _ -(l-|-2)Aobs .g, 

df ^ fiicmHiz) ^ H{z) ’ ^ 

•where D{z) is the comoving angular diameter distance, and 
H{z) is the Hubble rate. A foreground signal that exists in 
{dx, 9y, f)-space can be characterized by considering its power 
spectrum, which lives in the Fourier-conjugate space (u, v, t), i.e. 
the (u, V, t) Fourier mode contributes an intensity I (Ox, dy, f) oc 
g27rj(ii0j,+t)e„+Ti/) ^ caveat, r is not exactly a local Fourier con¬ 

jugate of / (often referred to as rj) in a wide-field instrument except 
at the center of the field because u and v are dependent on / except 
at the center. Given the mapping of Eqs. l |7|8| l, a Fourier mode in 
the observable cube maps into a Fourier mode in the cosmological 
cube according to 


kx — 


2ttu 

W)' 


2'kv 


and kz 


—2nTH{z) 
(1 “b 2:)Aobs 


In the context of a radio interferometer, the angular conjugate 
space (u, v) corresponds to transverse baseline separation in units 
of Aobs, and r to time delay. (The identification of conjugate- 
frequency with time delay is a consequence of the fact that the 
power spectrum of the electric field from a source and its corre¬ 
lation function are Fourier transforms; see e.g. [Parsons et al.|2012| 
for a rigorous treatment.) A foreground that is spectrally smooth 
in / contaminates only the small (in absolute value) conjugate- 
frequencies T, and hence small kz = fey in the cosmological box. 

Since there are not many Fourier modes available for small fe||, 
removing the small-fc|| modes has only a small impact on the pre¬ 
dicted precisions. But while the main astrophysical foregrounds are 
indeed spectrally smooth, if an antenna has a far side-lobe response 
then the cross-correlation of the electric fields measures with two 
such antennas can have a contribution at long time delay r (or in 
frequency space, it can vary rapidly as a function of /). We summa¬ 
rize the derivation of the effect here, drawing extensively on deriva- 


tions in the literature JDatta et al.||2010 

IVedantham et al.||2012[ 

[Morales et al.|2012||Parsons et al.[2012l[' 

rott et al.|2012||Liu et al.[ 


|2014a|bt . 


For two antennas separated by distance L (assumed perpen¬ 
dicular to the target field, so that = L/Aobs), the cos¬ 

mological transverse wave number being sampled is 



2TVy/ V? + v'^ 


2nL 

Aobs 77 ( 2 ) 


( 10 ) 


A foreground in a far side-lobe produces a time-delayed signal 
in the two antennas according to the geometric propagation time, 
which may be as large as L/c, so a contaminated region in the 
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The foreground wedge and 21 cm BAO surveys 5 


Fourier domain is produced at |t| ^ LJc. Mapped from the ob¬ 
servable cube into the cosmological cube, this corresponds to 


= Ifc^l = 


2-k\t\H{z) ^ 2-kLH{z) _ D{z)H{z) 


(1 + ^)Aobs c(l -t- 2:)Aobs c(l -t- 2:) 


k^. 

( 11 ) 


Equation fuf defines a “wedge” in Fourier space that is con¬ 
taminated by the response to sources far from the main beam due 
to geometric time delay, and a conicaQ “window” in which this 
process produces no contamination. 

Some words of warning are in order before we proceed. The 
wedge does represent information that is fundamentally irre¬ 
trievable, since with excellent knowledge of the side-lobe response 
and maps of the radio sky one could compute the contribution from 
sources far from the main beam and subtract it off (or solve for the 
contaminating sources as part of a global model). The “windows” 
are also not perfectly clean: other effects, notably reflections in ca¬ 
bles or antennas, can lead to a sky signal received at time t being 
detected again at some later time t -f At, and thus producing spu¬ 
rious power at r = Af. A similar effect can be produced by linear 
polarization-to-intensity (Stokes Q,U —> 7) leakage in the pres¬ 
ence of Galactic Faraday rotation: the arrival times of the left- and 
right-circularly polarized pulse from a distant electron are sepa¬ 
rated by a relative delay Afb-R because of the different indices of 
refraction for the two circular polarizations, which produces con¬ 
tamination in the Fourier modes corresponding to 77 = Afn-R 


(e.g., |Haverkorn et al.||2003a|b| |de Bruyn et al.||2006| |Pen et al.| 

[200^ IBemardi et al.||2009| |2010[ [Moore et al.||201^ (c.f., for a 

removal technique, |Shaw et al.|2014|[2015^ . 

However, multiple recent 21 cm experiments, including the 
PAPER UPober et al.|[20T^ and MWA pillon et al.|M4i|Thya-| 
Igarajan et al. 2015^ instruments that have produced the tightest 
upper limit on the reionization signal, have found that the wedge 
effect is the brightest foreground-instrument interaction signal in 
their data, and is orders of magnitude brighter than thermal noise. 
Therefore in this paper we forecast the effect of the wedge on BAO 
measurements using the 21 cm intensity mapping technique. The 
basic question is to determine how much loss of precision occurs 
if the wedge is removed, or alternatively what gains would be real¬ 
ized with calibration techniques that clean out the wedge. In what 
follows, we address this question by considering BAO forecasts 
with two different cases - a “Default” case, assuming the full 3D 
Fourier domain is available and noise-limited, and a “Wedge” case, 
in which information in the wedge is assumed to be inaccessible. 


2.4 Incorporating the wedge effect into the Fisher matrix: 
case of no reconstruction 

It is straightforward to include the wedge in the BAO Fisher matrix 
in the absence of reconstruction. We simply incorporate the effect 
of the wedge by excluding all the signal-to-noise contributions (Eq. 
[T]( from the Fourier modes at /r < /imin when calculating the Fisher 
matrix: 

_ fell _ D{z)H{z)/[c{l + z)] 

/^min — ! -- — ~ s 's 




fl + {D{z)H{z)/lc{l + z)]r' 


^ In 3D Fourier space {kx,ky,kz), Eq. 0 defines two cones, one at 
fc|| > 0 and one at fc|| < 0. Since power spectra satisfy the mathematical 
identity P{k) = P{—k), one normally considers only positive fc||. 


Equivalently, we may replace the response (i.e., window) function 
in Eq. with 

iV(fc) = IFinst(fc)0(/r - /imin), (13) 

where /Tmin is the cosine of the wedge angle from Eq. In the 
absence of reconstruction, the two formalisms are equivalent, but 
when we study reconstruction, the window function formalism will 
be more useful due to its flexibility in possibly softening the wedge 
after reconstruction. For simplicity, we assume the same /imin cut¬ 
off before and after reconstruction. Again, note that the minimum 
cosine angle /imin is a purely function of cosmology and the red- 
shift of the aimed signal such as 77 (z) and D{z), not depending on 
the configurations of the telescope at all. Using our fiducial cos¬ 
mology, /imin = 0.55 at 2 ~ 1 and 0.76 at 2 ~ 2. Figure 
shows the /tmin as a function of redshift over our fiducial range 
0-7 < 2 < 2. Figureshows the effect of /tmin on the precision 
of the isotropic distance scale Dy as a sanity check assuming 
our fiducial redshift bin at 2 = 1 ; aov as a function of /tmin 
is roughly consistent with the expected dependence based on the 
number of Fourier modes that are lost, i.e., oc l/-yi — /tminl the 
deviation at high /tmin is probably because the effective damping 
length also increases as 'Ez becomes more important than Exy with 
increasing /t. Overall the wedge effect are consistent with the mode 
counting when other conditions of the survey is fixed. The blue 
points/line account for the triangular window function we defined, 
and the black points/line do not include any window function. The 
results using the window function from |Seo et aL| ( |2010b^ give a 
prediction in between the two cases. 

We include /imin ( 2 ) from Figure 1 in the Fisher forecast for 
our fiducial 21cm BAO survey and investigate the combined effect 
of the wedge for the survey. In Figure]^ the black and red squares 
show the expected precisions on Dv (i.e., ao^) without and with 
the wedge effect, respectively, based on the survey configuration in 
Tables 

When the wedge effect is included, we find that an^ increases 
about ~ 36 — 46 per cent over the fiducial redshift range. It is 
interesting that this ratio is almost constant over the entire redshift 
range even though the range of /imin varies from 0.446 (at 2 = 0.7) 
to 0.762 (at 2 = 2.0), which alone would imply a substantially 
more increase in aoy at 2 = 2 (i.e., based on fl — /imin, o-o^ 
is expected to increase 30 per cent more at 2 = 2 compared to at 
2 = 0.7). Moreover, E^/Sxy is proportional to the redshift-space 
distortions (e.g., |Eisenstein et al.|2007^ which is larger at high 2 ; 
a larger /tmin at high 2 would have to use a more damped BAO 
feature, thereby increasing ao^ ■ 

The main leverage comes from the redshift-space distortions 
that amplify the power of the Fourier modes along the line of sight 
more at high 2 ; for example, /3 at 2 = 2 is 0.94 while /3 at 2 = 0.7 
is ~ 0.47 assuming the clustering bias of HI is I. The discontinuity 
at 2 ~ 1.22 in Figure 3 with the wedge effect (red squares) is be¬ 
cause the fiducial telescope configuration changes near this redshift 
(see Table[TJ in a way that affects the modes along the line of sight 
more. 

The wedge effect we have investigated so far is a combination 


® Dy is calculated by deriving Fisher matrix of D{z) and H(z) and then 
contracting the Fisher matrix of D(z) and H{z) with the vector (1, —1) 
assuming radial and transverse distances scale the same way, as explained 
in |Seo & EisensteTn]j2007t . This Dy constraint is therefore equivalent to 
varying Dy with fixed Alcock-Pazynski jAlcock & Paczynski|1979t shape 
DaH, corresponding to a BAO distance scale derived from a spherically 
averaged clustering information. 
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Figure 3. Left: eiTor on Dy without (black) and with (red) the wedge effect for the fiducial survey. No BAO reconstruction is assumed here. Therefore, the 
black squares represent the results similar to those in|Seo et al.lpOlOb}. Right: nPo .2 of the fiducial survey. 
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Figure 1. The minimum cosine line-of-sight angle /Tmin sot by the wedge 
effect as a function of redshift given our fiducial cosmology (black squares); 
red crosses show the minimum angle assuming the Planck ACDM best 
fit cosmology ^Planck Collaboration et al.|2015^ . Note that wedge effect 
solely depends on the cosmological parameters rather than the telescope 
configuration. Assuming WMAPl, /imin = 0.55 at z = 1 and 0.76 at 
z = 2. 

of Fourier mode counting over the redshift bins of a reasonable 
fiducial survey without assuming the BAO reconstmction. We next 
include the BAO reconstruction in the presence of the wedge effect 
through the BAO damping scales. 


Figure 2. The effect of /Tmin on an isotropic distance scale Dy (squares) 
in comparison to what is expected based on the number of Fourier modes 
cut off as a function of /Tmin (lines). We assume a condition similar to 
our fiducial redshift bin at z = 1 except for the volume, i.e., we assume 
Exy = 5.93, = 10.8, P = 0.82, /cmax = 0.6/1 Mpc"! (in Eq.[^, 

andnPo .2 = 0.46 while assuming a volume of l/i~^ Gpc^. The reference 
line is the expected functional dependence of based on the number of 
Fourier modes cut off by /Tmin ■ Blue points/line include the triangular win¬ 
dow function effect, while the overlaid black points/line assume no window 
function. 

3 THE WEDGE EEFECT ON THE BAO 
RECONSTRUCTION 

3.1 BAO reconstruction in Eisher matrix calculations 


The BAO reconstruction technique was first introduced in |Eisen-| 
[stein et aT] j2007b[l and has been shown to be successful in sim- 


ulations as well as in observations (e.g., Seo et al.j 2008 jSeo 

et al.||2010a| Padmanabhan et al.||2009||Noh et al.||2009||Mehta 

et al. 

201 lljPadmanabhan et al.|2012l| Anderson et al.|2012l| White 

2015 

1 . The method assumes a linear continuity equation to estimate 
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the large-scale displacement field from the gravitational field con¬ 
structed from the observed nonlinear density fields of mass tracers 
(Eq. [g. In Fourier space, the continuity equation becomes 

q{k) = -±S(k)&(k), (14) 

where i5(fe) is the Fourier coefficient of the nonlinear tracer density 
field and q(fe) is the Fourier coefficient of the 3-dimensional dis¬ 
placement field. S(fc) is a smoothing kernel that we convolve with 
the observed density field to reduce the small-scale nonlinearities 
before imposing the continuity equation. 

Displacing the tracers back to the estimated original positions 
using the derived displacement field, which are then followed by 
additional simple steps, substantially removes the nonlinear mode¬ 
coupling effect on BAO (especially, mode-coupling between large 
and small wave modes) and reduces the damping of the BAO. We 
refer interested readers to |Eisenstein et al.] ( |2007b) and |Padmanab-| 
|han et al1 ( |2009| l for more details. 

The Fisher matrix calculation can include the improvement 
due to the BAO reconstruction technique through the nonlinear 
damping scale, Sxy and E^, in Eq. 0. Depending on the details 
of the reconstruction, the level of power spectrum amplitude rel¬ 
ative to the shot noise may change, but we ignore such effect in 
this study. We follow the exactly the same procedure that is sum¬ 
marized in |Mehta et al.|p011^ , which was designed to recover the 
linear redshift-space distortions after reconstmction rather than re¬ 
moving it. 

In galaxy surveys, the conventional assumption is that the 
BAO reconstruction halves Exy and E^ (e.g., |Seo et al.|2008) when 
the signal per Fourier mode is much larger than the shot noise, i.e., 
nP 3> 1 at fe ~ 0.2h Mpc“^ which is the scale relevant to the 
BAO feature. We define nPo .2 in an intensity-mapping radio sur¬ 
vey to be the ratio of HI clustering power to the sum of both shot 
and instrument-thermal noise at this scale: 

p ^ Pmjk = 0.2h Mpc~\fj. = 0) 

Pth(fc = 0 . 2/1 = 0)-I-Pshot 

A low signal-to-noise per mode nPo .2 ^ 1 means that the den¬ 
sity field is noisy and therefore the displacement field derived from 
Eq.tg would be noisy as well, thereby decreasing the efficiency 
of the BAO reconstruction. Note that we continue to use this defi¬ 
nition in the presence of the wedge, even though the fj, — 0 modes 
are in the wedge and hence excised. 

The right panel of Figure shows nPo. 2 , i.e., the signal to 
noise per mode at fc = 0.2h Mpc“^ for the fiducial survey as¬ 
sumed. The figure shows that nPo .2 < 1 over the entire redshift 
range. [Seo et al.| ( |2010b^ made a conservative decision to assume 
that BAO reconstruction is not efficient when a redshift bin has a 
signal to noise ratio less than 2. We do not make such assumption 
in this paper, but instead directly estimate the reconstruction effi¬ 
ciency at this low signal-to-noise limit using simulations. 


3.2 BAO reconstruction of the wedged field 

The wedge effect, when not corrected for, would substantially re¬ 
move the information on the underlying density field, thereby de¬ 
grading an estimation of the underlying displacement field as well. 
Therefore, the effect of the wedge is not merely a cut-off of the 
contribution from the modes with jj. < fimin in Eq. ( |13| l, but should 
also include the degraded efficiency of the BAO reconstruction by 
using different assumptions of Exy and Ez as a function of /j-min 
after BAO reconstruction. 


While the damping scales before the BAO reconstruction are 
straightforward to calculate analytically using the expected linear 
theory power spectrum, it becomes nontrivial to predict the corre¬ 
sponding values analytically as a function of noise level after re¬ 
construction, especially in the presence of the wedge effect. |Pad-| 
[manabhan et al.| ( [2009^ has explained the BAO reconstruction tech¬ 
nique using the Lagrangian Perturbation Theory (hereafter “LPT”) 
formalism. Using the formalism, in principle, an analytic approxi¬ 
mation of the damping behavior after the BAO reconstruction could 
be derived. [White| ( [^ 1 0| l applies the formalism to estimate damp¬ 
ing scales in the presence of various galaxy shot noise levels; how¬ 
ever, these analytic calculations have not been confirmed in detail 
against the damping scales from simulations and also no redshift- 
space distortions were considered. Instead of including the wedge 
effect in the analytic formalism, in this paper, we directly estimate 
the values using A^-body simulations. 

Note that the Gaussian damping function in Eq. 0 is an ap¬ 
proximation for the linear information that survived in the late-time 
density field, which is often referred to as a “propagator” defined in 
Eq.([g:i.e., a cross-correlation between the late-time density field 
and the linear density field relative to the linear input power spec¬ 
trum. We therefore can estimate Exy and Ez by directly estimating 
the propagators as a function of redshift as well as a function of the 
noise level. 


(Sun(k,z)SMn(k',z)) 


( 16 ) 


where Sun(k,z) is the initial linear fields that is linearly scaled 
to 2 , and 5{k,z) is from the final density fields at 2 . Then the 
BAO signal in the nonlinear power spectrum can be modeled as 
C{k,^,z)^Piin{k,fi,z). 


3.3 N -body simulations to calculate propagators 

Since the propagator depends on a correlation between the initial 
field and the final field generated with the same random seed, the 
sample variance associated with the calculation can be substan¬ 
tially reduced with only one realizatioip] We use one realization of 
Particle-Mesh dark matter simulations that were used in ISeo et al.l 
( |2010a| >. The fiducial cosmological parameters were chosen to be 
consistent with WMAP5-I-SN-I-BAO results ( [Komatsu et al.|200^ : 
nm = 0.279, Da = 0.721, h = 0.701, = 0.0462, Us = 0.96, 

and as = 0.817. The simulation volume covers 2^h~^ Gpc^. The 
density field is derived using the Cloud-in-Cell (CIC) methodp^for 
576^ meshes that cover the volume of 2^h~^ Gpc® and Fourier- 
transformed. 

Note that this cosmology is not identical to the fiducial cos¬ 
mology used for the Fisher matrix calculations. We pick two red¬ 
shift outputs of one realization, at 2 = 1 and 2 = 2. We assume 
tracers with bias of unity, as shown in Table [T] by simply taking 
dark matter particles while varying levels of signal-to-noise by sub¬ 
sampling the dark matter particles. Although it is possible that, at 
a fixed bias, propagators for tracers with different HODs can be 


® We also apply Savitzky-Golay smoothing when plotting the propagators 
to reduce the noise due to the thin fc-shell averaging. 

In the Cloud-in-Cell method, we decrease the discreteness effect of a 
mass tracer particle by treating each particle to be a uniform-density cubical 
cloud with a volume of (2000/576)^h“® Mpc® centered at the particle 
position and count the mass fractions contributed to nearby 8 meshes. 
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slightly different from each other, the difference has shown to be 
small and therefore we ig nore it (e.g.|Padmanabhan & White|2009[ 
[Padmanabhan et al.|2009[[Mehta et al.|201 1^. 

To simulate different values of nPo. 2 , we subsample the 
dark matter particles at z = 1 and 2 in the following way. At 
z = 1, Po. 2 , i-e., the power at fc = 0.2h Mpc“^ is expected 
to be Mpc® using our fiducial cosmology. This requires 

the particle number density of ~ 0.007, 0.0028, 0.0014, 0.0007, 
and 0.00014/1® Mpc“® to mimic nPo .2 = 5, 2, 1, 0.5, and 
0.1 respectively. The simulation has a particle number density of 
0.0239/1® Mpc“®, and therefore we subsample particles by a frac¬ 
tion of 0.29, 0.117, 0.0586, 0.029, 0.00586, respectively. Note that 
we assume that this Poisson subsampling noise represents the noise 
due to the instrumental and sky background as well as shot noise, 
Pth{^') P T^shoti itt Eq. 0- We do not explicitly include the effect 
of the response function to represent a suppression of signal rela¬ 
tive to the noise on small scales. The instrument beam is smaller 
or comparable to the size of the mesh we use for the density calcu¬ 
lation and the displacement field calculation, as well as the typical 
nonlinear scale. Since we expect little information from below the 
mesh scale, we expect this procedure to accurately capture the dis¬ 
tribution of signal-to-noise ratio as a function of k. 

At z = 2, linear Po .2 is expected to be 342/i“® Mpc®, which 
requires the particle number density to be 0.015,0.0058, and 0.003, 
0.0015, 0.0003 /i® Mpc"® to make nPo .2 = 5, 2, 1, 0.5, 0.1, re¬ 
spectively. We subsample particles in the simulations by a fraction 
of 0.628, 0.243, 0.126, 0.0628, and 0.0126, respectively. 


3.3.1 Implementing the wedge 

The minimum p, angle that defines the wedge effect is approxi¬ 
mately 0.55 at 2 = 1 and 0.76 at z = 2 using our fiducial cosmol¬ 
ogy of the Fisher matrix analysis and we apply the same value for 
the simulation despite the small difference in cosmology. To mimic 
the wedge effect, we assume that Fourier modes for p < pmin are 
completely dominated by noise and simply set the corresponding 
Fourier coefficients to be zero. 

The pre-reconstructed power spectrum in the presence of the 
wedge is straightforward. The squared amplitudes of the Fourier 
coefficients of the pre-reconstructed density field (after corrected 
for the CIC pixel-window function) is recorded as a function of p. 

In order to derive a post-reconstructed power spectrum with 
the wedge effect, we have to first simulate a pre-reconstructed den¬ 
sity field in the presence of the wedge effect, hereafter and the 
displacement field that is derived from This is done in three 
stages: (1) the raw pre-reconstructed density field is Fourier trans¬ 
formed; (2) the wedge is applied in Fourier space by zeroing all the 
resulting Fourier coefficients when p < pmin', (3) we then inverse- 
Fourier transform these wedged Fourier coefficients to construct 
po,w. Note that the resulting density field is not particle-based any 
longer, but mesh-based. The displacement field with the wedge ef¬ 
fect is derived in a very similar way, but by modifying the second 
and the third stages: i.e., (2) by deriving the wedged q{k) (Eq. |14| 
by zeroing all q(k) when p < pmin and (3) then by an inverse 
Fourier transformation of q{k). 

We displace the mesh-based po,™ field according to the de¬ 
rived displacement fields; note that this is different from the con¬ 
ventional way of displacing galaxies, i.e., a particle field. We find 
that displacing a mesh-based density field is almost identical to dis¬ 
placing particles as long as the CIC pixel-window function is prop¬ 
erly corrected for the mesh-based density field. We then displace a 
uniform density field with the same displacement field and subtract 


the two in order to derive the final, reconstructed density field pr,™ 
in the presence of the wedge effect. 


3.3.2 Smoothing scale 


The BAO reconstruction depends on the smoothing kernel S{k) 
that we apply to the nonlinear density field to filter out high k 
modes when deriving displacement fields in Eq. 0- The opti¬ 
mal smoothing kernel for BAO reconstruction is shown by |Tassev| 
|& Zaldarriaga| ( |2012^ . However, the goal of the current particles- 
per is not to find an optimal BAO reconstruction, but to find an 
first-order description of the qualitative effect of the wedge, and 
therefore we assume a simple Gaussian smoothing kernel that has 
been conventionally used for simulations as well as for data. The 
choice of the smoothing length depends on the level of nonlinearity 
and also on the level of shot noise. [WhiteH2010l l shows that using 
a very small smoothing length when the high k is noise-dominated 
makes reconstruction ineffective while there exists a broad range of 
optimal/stable smoothing length. We want to choose the smooth¬ 
ing length that falls in this stably optimal range. We choose the 
smoothing scale at the noiseless limit to be 5/i“® Mpc (i.e., at ~ 
zero shot-noise limit) and this corresponds to the s moothing scal e 
of 5s/2h~^ Mpc using the convention defined in 


White 


1 2010 1 . 

We would want to alter this scale based on the shot noise level of 
subsampled field as shown Eq. l |17^ : 

1 


S(fe) 




1 -I- l/[nP{k )]' 


(17) 


where 1/(1 + l/[nP{k)]) is an optimal weight to filter out high- 
noise Fourier modes and derive the true S{k)5{k). Again, we 
choose = 5h ^ Mpc. 

For simplicity, however, we approximate this equation by 
defining an effective smoothing scale such that 


S{k) = e~ 


k^S%/2 

? 


(18) 


where E/j is set to produce 

_ g“^0.2^K,o/2 


1 

1 -f l/[nP{ko. 2 )] 


(19) 


at fc = 0.2h Mpc"^. This gives Er; = v/iM, v/45, v/OT, v/^, 
and \/145 h~^ Mpc, for nPo .2 = 5, 2, and 1, 0.5, and 0.1, re¬ 
spectively. Figure]^ shows the difference between Eq.[T^(red) and 
[T^(blue) which quickly becomes worse for low nPo. 2 . We tolerate 
this difference. 
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Figure 4. The difference in the smoothing kernel model between Eq. [m 
(red) that incorporates the optimal weighting and Eq. |18| (blue) that we ac¬ 
tually use as approximations. 

4 RESULTS 

In this section, we estimate the propagators at z — 1 and 2 for 
various values of nPo .2 using the A^-body simulation and adopt the 
derived propagators to model the BAO signal in the Fisher matrix 
predictions for 21cm radio BAO surveys. 

4.1 Reconstruction at 2 = 1 

Fi gure [5] shows the propagators before (black lines for nPo .2 = 
5) and after reconstruction (blue, magenta, and red lines for 
nPo .2 = 5, 1, 0.1, respectively), without (top) and with this wedge 
effect (bottom panels). The left panels show the propagators for the 
transverse modes (fj, — 0.05 ± 0.05), the middle panels show ones 
for the modes that are just outside of the wedge (/r = 0.65 ± 0.05), 
and the right panels show ones for the modes along the line of 
sight (/r = 0.95 ± 0.05). Note that we use modes slightly away 
from /r = 0 and /r = 1 in order to increase the number of 
modes and therefore to decrease sample variance in estimating 
Exy and S 2 from the propagators. The difference in jj,^ caused by 
An = 0.05 is small, but is accounted for when we derive Exy and 
E 2 . When estimating Exy and E^ for a Gaussian damping model, 
we use the propagator values at k = 0.3h Mpc“^ rather than at 
k = 0.2h Mpc“^ because, based on our visual inspection, the for¬ 
mer seems to give a better description over the scales of our interest 
(i.e., k < 0 . 3/1 Mpc“^) than the latter and the larger number of 
modes at the bin centered at k = O.Sh Mpc“^ allows the estima¬ 
tion of the damping length to be more stable. To further decrease 
the noisy fluctuations in the propagators, we smooth the measured 
propagators using the Savitzky-Golay filtering. The colored dotted 
lines correspond to the estimated damping length Eni,eff at given n 
based on the propagator value atk = 0.3h Mpc“^. 

In the presence of the wedge (bottom panels), the improve¬ 
ment in the modes outside of the wedge (bottom middle and bot¬ 
tom right panels) by the BAO reconstruction (the difference be¬ 
tween the colored lines and black lines) is clearly visible, while it 
is less than in the wedge-less case (the corresponding top panels). 
The modes in different directions are coupled and the divergence 
operation is non-local; nevertheless, we find that the linear den¬ 
sity field fluctuations along the line of sight is still reconstructable 

Without BAO reconstruction, propagators for different values of nPo .2 
are very similar and the difference is dominated by the noise. 


to some extent even in the near absence of the transverse infor¬ 
mation. For example, for nPo .2 = 5, we derive Exy and E^ to 
be 5.16h~^ Mpc and 10.16/i“^ Mpc, respectively, before recon¬ 
struction and 1 . 92 / 1 “^ Mpc and 4.86/i~^ Mpc after reconstruc¬ 
tion without the wedge, by extrapolating from the measurements of 
the damping scale at ^ = 0.05 and 0.95 to fj, ~ 0 and 1. With the 
wedge effect, Exy and E 2 are extrapolated from the measurements 
of the damping scale at /r = 0.65 and 0.95; we find 6.12/i“^ Mpc 
for Exy and 5.84/i“^ Mpc for E^. The decrease in E^ is still sub¬ 
stantial even with the wedge effect. One could consider generating 
constrained Gaussian realizations to fill the missing information in¬ 
side the wedge in an attempt to alleviate the effect of the wedge 
on reconstructing the modes around the wedge ([Hoffman & Ribakj 
|1991||Padmanabhan et al.|2012[ l, while our current results show a 
rather conservative scenario. 

The figure shows that the reconstruction efficiency decreases 
rather slowly when the signal to noise decreases from nPo .2 = 
5 to 1, but the efficiency decreases more drastically when nPo .2 
decreases from 1 to 0.1. Note that nPo .2 = 1 corresponds to the 
highest signal-to-noise redshift bins in our fiducial survey near 2 ~ 
1 (Figurej^. 

The degradation of the reconstruction efficiency in the pres¬ 
ence of the wedge becomes less obvious as the signal-to-noise level 
nPo ,2 decreases to 0.1, i.e., the largest noise case we investigate. 
The effect of the wedge on the BAO reconstruction is small in this 
case. It is probably because the modes within the wedge are already 
dominated by the noise. 

4.2 Reconstruction at 2 = 2 

A density field at a different redshift would respond differently to 
the combination of reconsbuction and the wedge. We expect the 
wedge situation to become more serious at higher 2 since a larger 
fraction of the modes are lost. Figure [^ shows the propagators be¬ 
fore and after reconstruction with and without this wedge effect 
for the transverse modes (n = 0.05 ± 0.05), modes that are just 
outside of the wedge (fi — 0.85 ± 0.05, middle panels), and the 
modes along the line of sight (/r = 0.95 ± 0.05). As expected, the 
reconstruction efficiency is significanly reduced in the presence of 
the wide wedge with n ^ 0.76. Again Exy and Ez are derived 
using the measurements at fi = 0.05 and 0.95 without the wedge 
effect, but derived from jj. = 0.85 and 0.95 with the wedge effect. 
Due to the larger sample variance in the two reference fi bins of the 
latter case, the derived estimates of Exy and Ez would be more sus¬ 
ceptible to noise. We again find that non-zero improvement in the 
line-of-sight modes outside of the wedge by the BAO reconstruc¬ 
tion, but at the level much less than that of the wedge-less case; the 
degradation due to the wedge is more severe compared to the cases 
at 2 = 1 due to the line-of-sight modes being closer to the wedge 
boundary. Again, the degradation due to the wedge is less obvious 
as the noise level increases to nPo .2 = 0.1; however, the effect 
of the wedge appears more severe even at this noise level com¬ 
pared to the same case at 2 = 1, probably because more modes are 
removed. Note that nPo ,2 = 0.1 corresponds to the lowest signal- 
to-noise redshift bins in our fiducial survey populated near 2 ~ 2 
(Figure]^. 

4.3 Propagators as a function of signal-to-noise and the 
wedge angle 

Table [^ summarizes the nonlinear damping scales for a Gaussian 
approximation we derived, with and without the wedge effect, and 
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Figure 5. Propagators of the subsampled redshift-space matter field, i.e., nPo .2 = and 0.1, at 2 : = 1 of modes in the transverse direction (left), modes 
along the line of sight (right), and modes just outside of the wedge (middle). Top: without the wedge effect. Bottom: with the wedge effect with //min = 0.55. 
Solid lines: measured propagators. Dotted lines: Gaussian model propagators based on the estimated damping scale at given // based on A: = 0.3h Mpc~^ 
(i.e., corresponds to the quoted Sni,eff)- Black lines: before reconstruction using the density field with nPo .2 = 5. Blue, magenta, and red lines: after 
reconstruction for nPo .2 = 5, 1, and 0.1, respectively. The gray reference lines are at unity and at the linear redshift-space distortion prediction. 



0,001 


0.01 0.1 
k (h Mpc^*) 




Figure 6. Propagators of the subsampled redshift-space matter field, i.e., nPo .2 = 5, 1, and 0.1, at 2 : = 2 of modes in the transverse direction (left), modes 
along the line of sight (right), and modes just outside of the wedge (middle). Top: without the wedge effect. Bottom: with the wedge effect with //min = 0.76. 
Solid lines: measured propagators. Dotted lines: Gaussian model propagators based on the estimated damping scale at given // based on A: = 0.3h Mpc~^ 
(i.e., corresponds to the quoted Sni,eff)- Black lines: before reconstruction using the density field with nPo .2 = 5. Blue, magenta, and red lines: after 
reconstruction for nPo .2 = 5, 1, and 0.1, respectively. The gray reference lines are at unity and at the linear redshift-space distortion prediction. 


© 2011 RAS, MNRAS 000,[T]-?? 




























The foreground wedge and 21 cm BAO surveys 11 


Table 3. The fit parameters ao and ai in Eq. j20). 


Case 

^xy 

ao 

ai 

ao 

ai 

2 = 1, no wedge 

+0.737 

0.659 

+0.562 

0.319 

z = 1, wedge 

-0.217 

0.837 

+0.453 

0.290 

2 = 2, no wedge 

-1-0.792 

0.773 

+0.619 

0.309 

2 = 2, wedge 

-0.544 

0.207 

+0.312 

0.291 


with and without reconstruction at z = 1 and 2. The last column 
shows the damping scale without reconstruction; it was evaluated 
using a set with nPo .2 = 5 since the nonlinear degradation be¬ 
fore reconstruction in principle would not depend on the noise level 
nPo. 2 - As a sanity check, we analytically calculated the expected 
damping scale without reconstruction using the method in|Eisen-| 
|stein et aL]j2007a| l and we derived very similar values (within 5 per 
cent). 

Figure|^visualizes the results in Table|^ i.e., the improvement 
in the nonlinear damping scale due to reconstruction as a function 
of nPo .2 in the presence of the wedge effect (red) and no wedge ef¬ 
fect (black); the damping scales are shown as ratios with respect to 
the values without reconstruction, i.e., Exy,o and Ez,o, and there¬ 
fore a smaller value corresponds to a bigger improvement. We find 
that the dependence of the reconstruction efficiency on nPo .2 is 
moderate from nPo .2 = 5 to nPo .2 = 0.5 (also see [Mehta et al.| 
[2011| | while it is substantial when approaching nPo .2 = 0.1. That 
is, we do not seem to gain much in terms of reconstruction by im¬ 
proving nPo .2 > 2. 

From the left panel of Figure]^ the smaller difference between 
the red and the black squares visualizes our earlier finding that the 
reconstruction efficiency for the modes along the lines of sight is 
not greatly affected by the presence of the wedge at 2 = 1, while 
the large difference between the red and black circles shows that the 
efficiency degradation when moving away from the line of sight is 
very large. The right panel visualizes the result at 2 = 2, i.e., a 
much larger degradation in the reconstruction efficiency due to the 
presence of the larger wedge. The differences between the red and 
black points decrease as nPo .2 decreases to 0.1, which again shows 
that the effect of the wedge is less as the noise level increases. The 
solid lines are fits to the ratios of the measured nonlinear damping 
scales assuming a fitting function in a form of 


Si f nPo,2 

Si.o V 1 + nPo.2 


ai 


( 20 ) 


Here i = xy or z is the direction (transverse or radial) along 
which we consider nonlinear smearing. The parameter ao describes 
the fractional improvement in the nonlinear smearing in the limit 
of high signal-to-noise ratio (nPo .2 ^ 1), whereas ai describes 
how that improvement depends on nPo .2 (ai = 1 would mean 
that half of the improvement is realized at nPo .2 = 1, whereas 
ai < 1 means that more than half of the improvement is realized 
at nPo .2 = 1). The fit parameters are shown in Table|^ 

For the radial nonlinear smearing length Ez, the main effect 
of the wedge is a modest reduction in ao. That is, if the wedge is 
excised, the reconstruction algorithm performs with the same de¬ 
pendence on nPo .2 (ffli ~ 0.3) but the amount of improvement 
(ao) is reduced by 20 per cent at 2 = 1 and 50 per cent at 2 = 2. 

The transverse nonlinear smearing length Exy exhibits more 
complex behavior. When there is no wedge effect, the dependence 
of Sxy/Exy,o on nPo .2 is stronger than that of Sz/Ez.o (larger 
ai), such that Exy/Exy,o is smaller than Sz/Ez,o at high values 


of nPo. 2 . This makes sense because the redshift-space distortions 
increase the effective nP(k) for the more radial modes that lead to 
displacements in the 2 -direction. The fact that Exy/Exy,o is larger 
than unity for the wedged case (ao < 0; see also the red circles in 
Fig.[8j appears at first sight to suggest that the reconstruction has 
made the nonlinearity problem worse in the transverse direction. 
However, we should recall that for a general Fourier mode with di¬ 
rection cosine fi, the effective nonlinearity scale is neither the trans¬ 
verse scale Exy nor the radial scale Ez, but rather the combination 
VSiy(l-/r2) + Ei/r2. That is, Exy is merely a parametric pre¬ 
sentation of the damping scale for the modes in the “dirty wedge” 
after reconstruction, extrapolated to /i = 0. The actual damping 
scale in the “clean window” (large y) decreased after reconstruc¬ 
tion, as shown in Figures[^and|^ 


4.4 The effect of the foreground on small line-of-sight wave 
modes during the BAO reconstrnction 

We have mentioned that the first order impact of the spectrally 
smooth astrophysical foreground would be losing the very small 
line-of-sight Fourier modes due to its contamination, and thereby 
we focused on the second order impact such as the foreground 
wedge that is more extensive in Fourier space. Again, the removal 
of the small line-of-sight modes has little effect on the BAO Fisher 
error forecasts due to the small number of modes, and therefore 
we have ignored these fc|| cutoff and assumed that the fluctuation in 
foreground contributes the additional noise that corresponds to lOK 
across all Fourier modes. However, the statement above warrants 
some caution when we consider the BAO reconstruction. While 
there are not many modes on large scales, the differential mo¬ 
tions such as the displacement fields are substantially contributed 
by large-scale modes bel ow k < O.lh Mpc~ ^ and even below 
k < 0.04/i Mpc“^ (e-g-, Eisenstein et al. 


2007b I. In this subsec¬ 


tion, we test the effect of ignoring the low fey cutoff in deriving the 
displacement fields. 

[Chang et al.H200^ presents very conservative limits for the 
smooth foreground subtraction, in terms of the wave numbers 
where the residual of the foreground subtraction due to the uncer¬ 
tainty in the foreground modeling exceeds the cosmological sig¬ 
nal. Their estimate gives k ~ 0 .01/t~^ Mpc for z ~ 1 and 
k ~ 0.02/i“^ Mpc at 2 = 2. As| Seo et af X 2(H0b| argues, this 
is a quite pessimistic limit. [Liu & Tegmark[ ( [2012^ , although they 
focused on the Epoch of reionization, shows that the foreground is 
featureless enough that it can be characterized to a great accuracy 
using only three or four independent parameters. Assuming that 
we need four Fourier modes to determine four independent param¬ 
eters for the foreground at each transverse location, we can derive 
fcz.min = 4 X 7r/Ar, where Ar is the radial range of the survey: 
for N polynomials, we would need fcz,min = {N -|- 1) x yr/Ar. 

For the fiducial survey in this paper, Ar ~ 2000h~^ Mpc 
between 0.7 < 2 < 2.11, which is equivalent to 0.60 e-folds in 
frequency. We choose fez,min = O.Ol/i Mpc“^ that corresponds to 
six (= N + 1) free parameters. Since our simulation box has the 
same dimension as Ar, this corresponds to removing six line-of- 
sight Fourier modes in the simulation box as well. When we remove 
all the modes for fez,min ^ O.Ol/i Mpc”'^ at 2 = 1 from the den¬ 
sity field (and therefore from the displacement field) in the case of 
nPo .2 = 5, i.e., the most effective reconstruction case, we find that 
the line-of-sight damping after reconstruction, Ez, increases by 16 
per cent, such that Ez/Ez,o in Figure|7]increases from 0.57 to 0.66. 
Meanwhile, there was little impact on Exy, which makes sense 
since the wedge has already removed most of the information in the 
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Table 2. Wedge effect on the non-linear damping scale of BAO. In the case without the wedge effect, in order to reduce the sample variance effect, Sxy and S 2 
are extrapolated from the damping scales measured at 0 < /i ^ 0.1 (centered at = 0.05) and 0.9 < /i < 1 (centered at ^ = 0.95), rather than measuring 
damping scales directly from Fourier modes at ^ = 0 and 1. In the case with the wedge effect, Exy and S 2 are extrapolated from the damping scales measured 
at 0.6 < fi ^ 0.7 (centered at /i. = 0.65) and 0.9 < fi < 1 (centered at = 0.95) for z = 1, and from the damping scales at 0.8 < ^ 0.9 (centered at 

fi = 0.85) and 0.9 < < 1 (centered at = 0.95) at z = 2. Rather than using a formal fitting procedure, we model the propagator as a Gaussian function 

with Exy and E 2 evaluated at k = 0.3/i Mpc“^. The damping scale without reconstruction in the last column is evaluated using nPo .2 = 5. 






with reconstruction 




without reconstruction 


nPo .2 

Sxy 

= 0.1 
s. 

nPo .2 

Sxy 

= 0.5 

s. 

nPo.: 

Sxy 

3 = 1 

s. 

nPo.: 

Sxy 

3 = 2 

s. 

nPo.: 

Sxy 

3=5 

s. 

Sxy,0 

Sz,0 

2 — 1, flulin — 6 

4.69 

7.63 

3.09 

5.96 

2.61 

5.51 

2.28 

5.18 

1.92 

4.86 

5.16 

10.16 

2 = 1, /J-niin ~ 0.55 

5.35 

7.93 

5.55 

6.72 

5.81 

6.35 

5.97 

6.10 

6.12 

5.84 



2 = 2, /J-min — 6 

3.36 

5.16 

2.17 

3.99 

1.92 

3.65 

1.53 

3.29 

1.16 

3.04 

3.59 

7.26 

2 = 2, /imin = 0.76 

4.63 

6.17 

5.44 

5.56 

5.21 

5.43 

5.40 

5.24 

5.36 

5.14 






Figure 7. The effect of the wedge on reconstructed nonlinear damping scales from Table The damping scales are shown as ratios with respect to the 
values without reconstruction, i.e., Exy,o and E 2 ,o- A smaller ratio below unity indicates more efficient BAO reconstruction. The black points show the 
results without the wedge effect and the red points show the results with the wedge effect. Circles show Exy/Sxy,o. i-^., the damping scale ratio for the 
transverse modes and squares show Taz/Taz,o^ i-^., the damping scale ratio for the line-of-sight modes. The lines ai‘e our models for the data points, i.e., 
1 — 0.0(^^0.2 /(I + ^^0.2))“^ T and the best fit oq and oi in each case are described in the main text. 
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Figure 8. The effect of wedge on the isotropic distance precision as a function of nPo. 2 - We assume a volume of 1 h~^ Gpc^ rather than the actual 
volume of each redshift bin of the fiducial survey. We assume the transverse resolution limit for the fiducial survey applies here such that the maximum 
^max = 0.6/i Mpc~^ at 2 = 1 (squares) and 0.35 at 2 = 2 (circles). Left: with reconstruction. Right: without reconstruction. The red points show the results 
with the wedge effect and the black points show the results without the wedge effect. 
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transverse direction that would be removed by In the lowest 

signal-to-noise case, i.e., nPo ,2 = 0.1, we find a smaller difference: 
an increase of 7 per cent in With k^^min = 0.007/i Mpc“^, 
which corresponds to four free parameters, we find 11 per cent in¬ 
crease in T,z with nPo .2 = 5. With fez,min = 0.007fe Mpc“^, 
which corresponds to four free parameters, we find 11 per cent in¬ 
crease in Ez with nPo .2 = 5. 

At 2 = 2, removing fez,min O.Olfe Mpc“^ gives Ez = 
5.47fe“^ Mpc after reconstruction for nPo .2 = 5, which is only 
6 per cent larger than the value in Table This value will in¬ 
crease Ez/Ez,o for nPo .2 = 5 in Figure|^from 0.71 to 0.75. The 
difference in Ez remains at 7 per cent at nPo ,2 = 0.1. Taking 
fez,min ^ 0.02fe Mpc“^ would increase Ez after reconstruction 
(for nPo ,2 = 5) by 13 per cent. 

We therefore conclude that our assumption of ignoring the loss 
of the very small line-of-sight Fourier modes due to the smooth 
foreground component is fairly good approximation even when 
considering the BAO reconstruction. 

4.5 Effect of the wedge on the distance measurements 

We next incorporate the resulting propagator estimates into the 
Fisher matrix calculation and estimate distance precisions while 
assuming our fiducial telescope configuration but only at redshift 
bins located at a = 1 and 2. In detail, we assume nPo .2 = 
0.1, 0.5,1.0, 2, 5.0 and a survey volume of lh~^ Gpc® in Eq. ([^ 
while taking a fiducial power spectrum and the resolution limit 
fcmax for at the corresponding redshifts defined by our fidu¬ 
cial configuration. Nonlinear damping scales Exy and Ez are taken 
from Tablej^given the redshift and the noise level. 

From the right panel of Figure]^ i.e., without reconstruction, 
one can see that the distance precisions at « = 2 is slightly better 
than those at z — 1 without the wedge effect. If the effect of the 
response function at the two different redshifts (i.e., femax in W^) 
were ignored, the distance precision at z = 2 would have been no¬ 
ticeably better due to the smaller nonlinear damping effect. With 
fcmax = 0 . 35/1 Mpc“'^ that limits along the transverse direc¬ 
tion at 2 = 2, compared to 0.6h Mpc“^ at 2 = 1, the precision 
of Da (z) at 2 = 2 becomes substantially deteriorated, making the 
overall precision of Dv only slightly better at 2 = 2 with help from 
a larger value of p. The 1 — a error contours of Da{z) and 74 in 
the right panel of Figure[^(black ellipses) shows that the precision 
of Da (z) constraints is slightly worse at 2 = 2 while the precision 
of H (2) is much better at 2 = 2. 

With the wedge effect (red points in the same panel), the er¬ 
ror on Dv increases by a factor of 1.27-1.8 (i.e., larger at higher 
nPo. 2 ) at 2 = 1 and 1.36-2.1 at 2 = 2 reflecting the more severe 
wedge effect at 2 = 2. As a result, the error on Dy with the wedge 
effect are coincidently almost the same at 2 = 1 and 2 = 2. In the 
right panel of Figure|^(red error ellipses), one can see that Da{z) 
is virtually not constrained at the 2 = 2 bin. 

The left panel of Figure]^ shown in comparison to the right 
panel, shows that the reconstruction is more effective at 2 = 1 ex¬ 
cept at the very low signal-to-noise limit. Figure|^implies a slightly 
more effective reconstruction at 2 = 2 in terms of the damping 
scale decreasing, but decreasing damping scale brings a larger re¬ 
turn at the lower redshift when the field is more nonlinear to begin 
with. By comparing the red points in the two panels, we calcu¬ 
late the reconstruction efficiency in the presence of the wedge; the 
fractional decrease in error is 0.85-0.65 at 2 = 1 and 0.93-0.79 
at 2 = 2. Therefore, despite the wedge effect, the survey would 
benefit substantially from the BAO reconstruction. 


The left panel of Figure]^ shows the 1 — a ellipses of Da{z) 
and H ( 2 ) with reconstruction. Note that the ellipses with the wedge 
effect (red) do not touch ellipses without the wedge effect (black) 
after reconstruction, which is distinguished from the cases (the 
right panel) without the reconstruction. This is reasonable in that, 
with reconstruction, the wedge effect does more than just wiping 
out the first-hand transverse information; in the presence of the 
wedge, the reconstruction of modes along the lines of sight is also 
affected by the missing transverse modes because of the worse fi¬ 
delity of displacement fields. On the other hand, we find that if we 
do not remove the transverse modes until after we derive displace¬ 
ment fields (i.e., when the displacement fields are correctly derived 
and the wedge is then applied), the ellipses of the two cases (red and 
black in the left panel) tend to touch each other at least at 2 = 1. 

Finally, we apply the propagator estimates to all the redshift 
bins of our fiducial survey which have different fiducial power 
spectrum amplitudes, resolution limits, nPo. 2 , and the survey vol¬ 
umes at different redshift bins that are included. The fiducial 21 cm 
BAO survey we assumed has nPo .2 ~ 1 at 2 ~ 0.7 and decreases 
to ~ 0.1 at 2 ~ 2 (Figure]^. However, Figure [7] shows a non- 
negligible level of reconstruction efficiency even at nPo .2 ~ 1- 
We calculate the nonlinear damping scales for the WMAPl cos¬ 
mology by taking Exy,o = 9.38fe“^ Mpcj^at 2 = 0, rescale 
Sxy,o(-2) with the growth function, and derive Ez accounting for 
redshift-space distortions. The reconstruction is incorporated by 
making a use of the fitting models for Exy/Exy,o and Ez/Ez,o, 
i.e., 1 — ao(nPo. 2 /(l + nPo. 2 ))“^ with the estimates of oq and ai 
from § |4.3[ Since we have propagator estimates for only two red¬ 
shifts, i.e., at 2 = 1 and 2, we apply ao and ai derived at 2 = 1 to 
the redshift bins between 0.7 and 1.5, and ao and ai from 2 = 2 to 
the redshift bins between 1.5 and 2. 

Focusing on 77(2) constraints (the right panel of Figure [I^, 
we find that, the reconstruction improved the precisions from the 
wedged field (compare the red solid and red open points) by 40 per 
cent near 2 ~ 1 and 10 per cent at 2 ~ 2. For the readers who 
are interested in comparisons to |Seo et al.|p010b^ , near 2 ~ 1, 
the reconstructed field with the wedge effect (solid red) would give 
a better precision on 77 ( 2 ) than the un-reconstructed field without 
the wedge effect (open black, i.e., similar to the default from |Seo] 
|et al.| ( [2010b^ except for the assumption of the response function), 
while the reconstruction helps little near 2 ~ 2 due to the very 
low nPo .2 of the reference survey. For Da(z) constraints (middle), 
as expected, the wedge effect is detrimental; nevertheless, the re¬ 
construction makes ~ 20% improvement in errors on Da{z) near 
2 ~ 1 despite the wedge effect. The left panel shows the resulting 
constraints on Dv- 

In order to estimate the overall effect of the wedge when BAO 
reconstruction is applied, one would want to compare the difference 
between the solid black points and the solid red points; the wedge 
effect increases the error on H{z) by a factor of 1.5-1.6, the error 
on Da{z) by 3-4.4 times. In addition, the wedge effect substantially 
increases the cross-correlation coefficient from the original 0.41 
(low 2)-0.44 (high 2 ) to 0.62-0.74, reflecting the increased degen¬ 
eracy due the poorly constrained Da{z). Reconstruction slightly 
increases the coefficient for the wedged case, probably due to a 
larger offset between the error on Da{z) and H{z) after recon¬ 
struction. 


The value is rather large because of a large erg (= 0.9) for this cosmol¬ 
ogy- 
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V°* 

Figure 9. The effect of wedge on angular diameter distances and bubble parameters as a function of nPo ,2 ■ We assume a volume of 1 h~^ Gpc® rather 
than the actual volume of each redshift bin of the fiducial survey. We assume the resolution limit for the fiducial survey applies here such that the maximum 
fcmax = OSh Mpc“^ at 2 = 1 (squares) and 0.35 h Mpc“^ at 2 = 2 (circles). Left: with reconstruction. Right: without reconstruction. The red ellipses 
show the results with the wedge effect and the black ellipses show the results without the wedge effect. One sees that the wedge effect is more severe if we 
consider BAO reconstruction. 
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Figure 10. The effect of the wedge on distance constraints predicted for the fiducial survey. The open symbols show distance constraints without the BAO 
reconstruction, and the sold symbols include the BAO reconstruction. Black points: default, without the wedge effect. Red points: with the wedge effect. 


5 CONCLUSIONS 

Foregrounds from Galactic and extra-Galactic origins in combina¬ 
tion with the chromatic response of the instruments are expected to 
contaminate a larger wedged volume in Fourier space around the 
transverse modes. In this paper, we investigated the effect of the 
contaminated wedge on the future 21cm radio BAO surveys con¬ 
sidering the BAO reconstruction technique. 

We first used A-body simulations to find the efficiency of the 
BAO reconstruction in the presence of the wedge and compared 
it to the wedge-less case for various noise levels. We found that 
the BAO reconstruction can still improve the BAO signal for the 
Fourier modes along the line of sight despite the near absence of the 
transverse information. The efficiency is sightly less in comparison 
to the wedge-less case, and the degradation increases quickly for 
modes away from the line of sight. Such degradation also increases 
as the redshift increases due to the larger wedge at higher redshift. 
The degradation becomes less obvious as the density field becomes 
noisier. 

We find that the dependence of the reconstruction efficiency 
on the signal to the noise level per mode, as parameterized with 
nPo. 2 , is moderate from nPo .2 = 5 to nPo .2 = 0.5 while it is 
substantially worse when approaching nPo .2 = 0.1. The ability to 
do reconstruction saturates at high nPo. 2 , with little improvement 
at nPo .2 > 2. The fiducial 21 cm BAO survey considered here 
operates in a low signal-to-noise regime with 1 > nPo .2 > 0.1 
at 1 < 2 < 2. The BAO reconstruction appears still valuable for 
this level of noise. We derive fitting formula for the BAO damping 
scales as a function of noise with and without the wedge effect. 


The overall effect of the wedge for our fiducial 21 cm radio 
survey was derived using the Fisher matrix analysis utilizing our 
prediction on the BAO reconstruction efficiency based on the find¬ 
ings written above. As expected, the effect of the wedge is detri¬ 
mental to the measurements of the angular diameter distance; the 
errors on angular diameter distances increased by 3-4.4 times. The 
wedge effect increases the error on H{z) by a factor of 1.5-1.65 
on the other hand. Again, the reconstruction was still effective even 
with the wedge. For H ( 2 ), we observed 10-40 per cent of a preci¬ 
sion improvement between 2 = 1 (nPo .2 = l,Aimin = 0.55) and 
2 {nPo .2 = 0.1, ^min = 0.76), which was up to 10 per cent less 
improvement compared to the wedge-less case. Even for Da{z), 
we still see 20 per cent improvement despite the wedge. 

In summary, we find that the calibration technique to lessen 
the wedge effect would be extremely valuable for obtaining com¬ 
pelling measurements of angular diameter distances from the 21 
cm BAO surveys with the advent of near-future BAO surveys in 
other wavelengths. For example, |Shaw et al.| ( |2014| |2015^ show 
that given the perfect knowledge of the telescope we can clean the 
foregrounds even with polarization well within the wedge down to 
fc|| ~ 0.02h~^ Mpc, which essentially would be the same as the 
wedgeless case. 
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